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Abstract. — Extending previous results we show that there are no Petrov type III 
space-times on which either the Weyl neutrino equation or Maxwell's equations satisfy 
Huygens' principle. We prove the result by using Maple's NPspinor package to con- 
vert the five-index necessary condition obtained by Alvarez and Wunsch to dyad form. 
The integrability conditions of the problem lead to a system of polynomial equations. 
We then apply Maple's grobner package to show that this system has no admissible 
solutions. 

RESUME. — En prolongeant des resultats precedents, on demontre qu'il n'existe aucun 
espace-temps de type III de Petrov sur lequel l'equation de neutrino de Weyl ou les 
equations de Maxwell satisfait au principe d'Huygens. Nous prouvons le resultat par 
utilisant le logiciel NPspinor de Maple pour transformer la nouvelle condition a cinq 
indices obtenue par Alvarez et Wunsch en composants de repere spinoriel. A partir des 
conditions d'integrabilite du probleme, on obtient un systeme d'equations polynomes. 
Nous employons done le logiciel grobner de Maple pour demontrer qu'il n'existe aucune 
solution admissible de ce systeme. 



1 Introduction 

This paper is the sixth in a series devoted to the solution of Hadamard's problem for 
the conformally invariant scalar wave equation, Weyl's neutrino equation and source-free 
Maxwell's equations. These equations can be written respectively as 

Uu + Iru = 0, (1) 
o 

V A V = 0, (2) 

V aA <Pab = 0- (3) 

The conventions and formalism used in this paper are those of Carminati and McLe- 
naghan jS] (CM in the sequel). All considerations in this paper are entirely local. Part of 
the calculations presented here were performed using the NPspinor package available in 
Maple [3IH1. 
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Huygens' principle is valid for eqs. (^Q), (j2J) and Q if and only if for every Cauchy 
initial value problem and every point xq in the 4-dimensional pseudo-Riemannian space 
V4, the solution depends only on the Cauchy data in an arbitrarily small neighborhood of 
S n C~(xo), where S denotes the initial surface and C~(xq) the past null conoid from xq 
^1[201E3]. Hadamard's problem for equations (Q), @ and ©, originally posed only for 
scalar equations, is that of determining all space-times for which Huygens' principle is valid 
for a particular equation. As a consequence of the conformal invariance of the validity of 
Huygens' principle, the determination may only be effected up to an arbitrary conformal 
transformation on the metric in V4 

9ab = e 2ip g a b , (4) 



where ip is an arbitrary real function. 

It is known that Huygens' principle is valid for (|T|). (J2J) and (j2J) on any space-time 
conformally related to the exact plane- wave |10[ 1131 121 j , with metric 



ds z = 2dv du + [D(v)z z + D{v)z z + e(v)zz]dv - Idzdz , (5) 

in a special coordinate system where D and e are arbitrary functions. These are the only 
known space-times on which Huygens' principle is valid for these equations. Furthermore, 
it has been shown I14[ I21j that these are the only conformally empty space-times on 
which Huygens' principle is valid. 

In the non-conformally empty case several results are known. In particular for Petrov 
type N, the following result has been proved [3J|IJ[2]: Every Petrov type N space-time on 
which the equations (OJ), HP and 0j satisfy Huygens' principle is conformally related to an 
exact plane wave space-time (5). For the case of a Petrov type D the following result has 
been established E2l Eh There exist no Petrov type D space-times on which equations 
0^ HP or (pJJ satisfy Huygens ' principle. 

In this paper we complete the analysis for Petrov type III spacetimes given by CM, 
in the case of equations (^Q) and (J2J). Our main result is contained in the following theorem: 



Theorem 1. There exist no Petrov type III space-times on which Weyl's equation 
|^ or Maxwell's equations J3]j satisfy Huygens' principle. 



2 Previous results 

Let ^ abcd denote the Weyl spinor. Petrov type III space-times are characterized by the 
existence of a spinor field o A satisfying 

^ABCDO C D = , ^ABCDO D + 0. (6) 

Such spinor field is called a repeated principal spinor of the Weyl spinor and is determined 
by the latter up to an arbitrary variable complex factor. Let l a be any spinor field satisfying 

o A i A = l. (7) 
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The ordered set {oa, ia}i called a dyad, defines a basis for the f-spinor fields on V4. The 
main results obtained by CM can then be stated as follows: 



Theorem 2. The validity of Huy gens' principle for the conformally invariant scalar 
wave equation (0), or Maxwell's equations J2J), or Weyl's neutrino equation J3Jj on any 
Petrov type III space-time implies that the space-time is conformally related to one in which 
every repeated principal spinor field oa of the Weyl spinor is recurrent, that is 

°A;BB = °Al B B . ( 8 ) 

where I BB is a 2-spinor, and 

^ABCD;EE^ B ^O D O E O^ = 0, (9) 

R = > ®abAb° A ° B = °- (10) 
Theorem 3. If any one of the following three conditions 

^ABCD;EE^ A ^ B ^^ E O E = 0, (H) 
*ABCD;EE<< A <' B O D O E l E = 0, (12) 
^ABCD;EE^ B ^ D O E O E = 0, (13) 

is satisfied, then there exist no Petrov type III space-times on which the conformally invari- 
ant scalar wave equation or Maxwell's equations or Weyl's neutrino equation J3]j 
satisfies Huygens ' principle. 

It will be proved here that conditions (fTT|) to (|13f) are superfluous in the cases of 
equations © and ©, i.e., they are consequences of the necessary conditions for the validity 
of Huygens' principle, in particular of the five-index necessary conditions derived by Alvarez 
and Wiinsch Q]|2]. 

3 Necessary conditions 

In order to prove the Theorem 1 we shall need the following necessary conditions for the 
validity of Huygens' principle P ITO ) ITfi | [TB I 120] : 

{III) S abk . k - \c k ah l L kl = 0, (14) 



(V) TS kiC k a b m Ckcdl t rn + ^k2C k a b-,cSkid + ^{8ki — k 2 )S a b k S c dk 

— 2k2C k abSklc;d — 8kiC k a b l S c dk:l + k 2 C k a b Q m c kLdm 



+ 4k\C k a b ''C m cdlLkm 



0, (15) 
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where 

Cabcd ■= Rabcd - 2flf[ a [d-Efft] c ] , (16) 

1 

Lab ■= -Rab + c R 9ab , (17) 
6 

S a bc '■= ^a[6;c] • (18) 

In the above C a \, c d denotes the Weyl tensor, R a b the Ricci tensor and TS[ ] the operator 
which takes the trace free symmetric part of the enclosed tensor. In (|15|) k\ takes values 3, 8, 
5 and k 2 the values 4, 13, 16, depending on whether the equation under consideration is the 
conformally invariant scalar equation, Weyl's equation or Maxwell's equations respectively. 

We shall also need a necessary condition involving five free indices, valid for Weyl's 
equation (2) and Maxwell's equations (3), that was found by Alvarez and Wiinsch ^HJ. It 
can be written in the form 

T W abcde + a x T^ ahcde + a 2 T^ abcde = , (19) 

where a\ and a 2 are fixed real numbers, and 

rp(l) r rQ\A*f~*k fi*/^^ 1 ^ {~ 1U _L f^ u k f^ v _1_ 

1 abcde — J O [4 U a b U del;ukc ~ O U be ;a ^ del;uk ~r ZD ^ ab ;u ^ dek;vc~r~ 

1 T 

cdl\k^en 



^ ao ; ^kdel\nc i" O O a fc U dei,c-^kn "T * l> ab ^ cdl;k^e 



where 



+ 4*C k a b C U Mc;u Lde ~ 2\* C k a b C U ' c dk;u L e l] , (20) 

th(2) T^CT 1 0*/"*& ^ n nfa */^ffc / 1 /r\-i \ 

J afcede •— J >J[— 1Z U a t, Ufc ca - U l ne - u ~ ^k aS^lnhb u de ;cj > (, Z1 J 

r7-i(3) T'OT q* r~ik I r~i n/-iu i * /^f n/i fik I l /oo\ 

1 abcde •— 1 o[—o O a fe Ufc c rf O £ ne;M + Cfc a^lbnh^ de ;c\ , 

Cabcd := 2 £a ' ,e ^efed- (23) 

It is worth to mention that, in our conventions, the Riemann tensor, Ricci tensor and 

the Ricci scalar have opposite sign to those used by Alvarez and Wiinsch 2 . The spinor 
equivalents of conditions III and V are given, respectively, by |171 @] 

*ABKL- K A L B + *AB KL $ KlA B = > ( 24 ) 

^1* ABCD;KK^ ABCD; KK + k 2 i & K (ABC;D)(A^ BCD)L; L 'K 
+ k 2 ^ K \abC;D)(A^ BCD)L- L k - 2 ( 8fe l - k ^ ABC\K[ K (A^ BCD)K- K D) 
~ k 2 ^ K (ABC * (ABC\L[ L '\K\D)D) ~ k ^ (ABC^ \ABC\L\? ' \K\D)D) 

+ Ak^ K {ABC * ' ABC]Ll . L D)Ki3) + ^ K ABd * ABClLl . L n )D)k 

+ 2(k 2 - 4k{)^ K {AB c^ D)K K(A^BCD) K ~ 2 ( 4A; 1 + h)^ABCD^ABCD = ■ ( 25 ) 
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4 Proof of Theorem 1 

In CM, Theorem 2 was proved by using conditions III and V. The explicit form of these 
necessary conditions is obtained by first converting the spinorial expressions to the dyad 
form and then contracting them with appropriate products of o A and i A and their complex 
conjugates. In particular, it was shown, that there exists a dyad {oa, la} an d a conformal 
transformation such that 

K = a = p = T = e = 0, (26) 

$0 = *1 = $2 = ^4 = 0, * 3 = -1 , (27) 

$00 = $01 = $02 = R = , (28) 

Da = Df3 = 0, (29) 

$n = c, (30) 

where c is a constant. By contracting condition (III) with i AB AB and l a o b l ab (where the 
notation 

OAi—Ap = OAi • • • oa p , etc. has been used) we get, respectively, 

L>vr = 0, (31) 

6/3 = -(3{a + (3) . (32) 
From the Bianchi identities, using the above conditions, we obtain 

D$ 12 = 2tT$ii , (33) 

£>$ 22 = -2(/3 + /3) + 2$2ivf + 2$i 2 7T, (34) 

5<S> 12 = 2a + 4vf + 2A$n - 2a$ i2 , (35) 

5$i2 = -2(3 + 27i$n - 2/3$i 2 , (36) 
5$ 22 = A$i2 + 27 + 4/7 - 2F$n + 2A$21 

+ 2$ 12 ^ - 2$ 22 /3 - 2$ 22 « + 27$i2 • (37) 
Using Ricci identities we get 

L> 7 = wa + (3tt + $11 , (38) 

5w = D A - ?f 2 - Wa + if (3 , (39) 

Dz7 = A7f + 7FTi+A7r + 7F7-7f7-l + $12 , (40) 

8a = 6(3 + aa + - 2/3a + $ u , (41) 

57r = Dp — Wtt + ira — (3tt . (42) 

We can obtain useful integrability conditions for the above identities by using NP commu- 
tation relations. Using (JHSJ, pEjl. (|37|) . (J32|) . (jSHJ, (EDI) and |@2J) in the commutator 
expression [5, Z?]$2 2 — [A,L>]<3?i2, gives 

5]3 = -2$n - /fo - 4/J7T - 2D/Z - # + 2vf7r. (43) 
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We shall consider both Maxwell and Weyl cases separately. We begin with the 
Maxwell equations, i.e., k\ = 5 and k% = 16. By contracting condition V with i / ABCD i AB o CD , 
we get 

62/?vr + 40/?a + 6vtq + 3a 2 + 5a + 25ir + 5/3 + /f = . (44) 
Substituting (|3*2*)) into this equation results in 

5(2W + a) = -62WP - 39(3a - 6Wa - 3a 2 . (45) 

Now, from (JUJ), (g^J and (@3J) we obtain 

<5(2vr + a) = 2vra + aa - 6/3vr - 30a - $ n . (46) 

Contracting condition V with i ABC o D t ABC o D , using (|38|). (|42[l and the complex conjugate 
of P3|). we get 

148$n + 152/?7f + 76(L>^ + Dp) - 8Wa - 104/?/? - 8aa - 232vrvr + 152/?vr - 8vra = . (47) 

Using (gBJ, (EU), gU, (El, (jH and gTJ) in [5,5](a + 2vr) = (a -^)<5(a + 2vr) + (-a + 
[3)5{a + 2tt), we obtain one expression for DJ2. Substituting it in l)47J) and solving for Z)^, 
we obtain 

Dfi = — f-1520a/3vr + 208aavr - 152/?7fvr + 1040/?/?a + 1968/?/3vr 
152 V 

+ 1216/?vra - 1228tt$h + 1688avr7f - 739a$n + 80a 2 a + 380/3$ u 
+24967T7T 2 - 2432/3tt 2 + 128vr 2 a + 80vfa 2 + 760/?aa) /(-/? + 5a + 8vr) , (48) 

where we have assumed that f3 — 5a — 8tt ^ 0. By substituting this equation into (|46|) . we 
find: 

DJI= — (776a/3a - 304/5 2 vf + 208/f /? + 84/J$n - 12167f7r 2 
152 V 

+ 741a$ii + 1536avf/5 - 760avfvr + 2744vfvr/3 + 13327r/?a 

+741a$n + 11407r$n) /(/? - 5a - 8vr) . (49) 

We note that (|48|) and ()49[) have the same denominator. So, in what follows we shall use 
the Pfaffians 5a, Sir, and 5(3, given by (|41jl. lj4*2*|) and (|43|) . respectively, and their complex 
conjugates, in such a way that they have all the same denominator. This procedure simplifies 
the expressions to be obtained from the integrability conditions. 

Now need to convert (|19[) to the spinor form in the dyad basis. The resulting 
expression, obtained using Maple's package NPspinor [H], has a considerable size, specially 
due to the term in (1201) containing a third order derivative of the Weyl tensor, and will not be 
presented here. However, among the twenty one independent spinor components, we found 
a relatively simple one, obtained by contracting the dyad expression with t/ AB P DE T ABGDE , 
It has the following form: 

- 14a 5w - 12a 3 - 5(5(a + 2w)) - 21/?a 2 + 7A/? 2 + U/3 2 W 

- 7/?<5(a + 27T) - 42/3a7f - (10a + 6Jf)5a - 24a 2 vf = . (50) 
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We observe that the terms (|21j) and (|22|) did not contribute to this component. Using (|45(1 
to eliminate 5tt from this equation, and solving for 5a we get 

5a = 192(3w - 3a 2 + 121/35 (51) 

and 

5W = - 127w{3 - 80/35 - 3Wa . (52) 

We have now all the Pfaffians we need to complete the proof. The integrability conditions 
provided by the NP commutation relations can now be used. Let us consider the NP 
commutator [5, 5]a. Using the Pfaffians calculated previously, and solving for <&n, we 
obtain 

$n := -8/3 ( 8 vr + 5 a ) (172736 tttt 2 - 7776 tt 2 a - 13294 air a + 7866 /3tt a 
+ 211556a7fvr + 9568^/3 tt - 22572 finn + 5330 /?^a + 4845 /3aa 

- 1805/3~ 2 a- 2470/3 /3~ 2 - 13110 na /3 - 5290 a 2 a + 65010 W a 2 )/ 

(-772320 tt 2 /?- 2048352 a tt 2 - 1158240 tt a 2 + 4085 of + 10868tt/3" 2 (53) 

- 335985 a 2 - 1020276 a tt /3 + 214700 a 3 + 1191680 tt 3 ) . (54) 

On the other hand, from the commutator \5, 6] (a + (3) the following expression for <E>n 
results: 

$n := -8/3(-20672/37fvr- 129207fa/3 + 8056/37ra + 1727367f7r 2 

- 7776 tt 2 a + 9568 /3 /3 tt + 5035 /3 a a - 5290 a 2 a + 65010 W a 2 

- 13294a7ra + 211556a7f7r + 5330^/3a)/(-13047^a (55) 
+ 42940a 2 + 162944 tt a + 10412 /3~ 2 - 18564 tt/? + 148960 tt 2 ) . (56) 

Using the fact that 5(3>ii) = 0, we obtain, from 1)55(1 . a third expression for 3>n: 

$n := 8/3(-48191081692160 7r 3 aa 2 + 180931104170496 tTtt 3 a 2 (57) 
+ 7968511840 /f 4 7t5- 265004094784 /f /3 tt 2 - 5268590200832 tttt 4 /3 (58) 
+ 24894675520 ~$ (3 tt + 221453789400 tt a 2 /f - 20739582848 ~f p 2 a 
+ 4980319900 ~f a a- 30330200378072 n 2 aa 3 - 1208861900450 a 5 a 
+ 52252921810507fa 5 + 724772046800 ~f tt tt a 

- 3675701240760 ^ 2 7f7r a 2 - 6145239989312^ tttt 2 a 

- 7649757648780 /3tt Tra 3 - 16323677160464 ]3tttt 2 a 2 

- 15274850502912 ]3tttt 3 a - 26601021440 /f tt a a 

- 1040399202440 /3 2 tt a a 2 - 1635701635136 /if tt 2 a a 
+ 5632482563850 /3"tt a a 3 + 13162176133400 /?tt 2 a a 2 
+ 13687783600768 /3tt 3 a a - 329130147840 /f /3tt a 
+ 265058659320 /f /3 tt a 2 + 469041329536 ~f (3 tt 2 a 

+ 1952913512680/3/3 na 3 + 4375379783424 /3/3 tt 2 a 2 

+ 4338917366784/3/3 tt 3 a - 102469061500 /f /3 a 2 + 49609754350 ~f (3a 3 

+ 325649974650/3/3 a 4 + 15559172200 /^ (3 a - 38360907652096 tt 4 a a 
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+ 39628423187260 7f a 4 n - 1329967209650 tt a 4 (5 - 858758431040 $ tt 3 a 

- 20447502080 tttt - 8524238850 /f a 2 a - 12779688800 Tajf 

+ 905003032125 /?a 4 a + 136081885849600 tttt 4 a 

+ 119915073751888 tttt 2 a 3 - 9564050953390 a 4 a tt 

+ 5344911334400 7r 4 5/3- 12237656133632 n 5 a + 1606267826176^/5 tt 4 

+ 40801870077952 tttt 5 - 730927683900 tt a 3 ~f - 220969445675 /f a 3 a 

+ 274949282816^ 2 /3tt 3 + 592713574016^ TfTT 2 -3416931669632 ^ 2 7f7r 3 ) / 

[(-13047/?a + 42940 a 2 + 162944 tt a + 10412 ]f - 18564 tt/3 + 148960 tt 2 )( 

102596352 tt 3 - 7684576 tt 2 ]3 + 194465152 a tt 2 + 24320 ir]f - 9460852 a ir~(3 

+ 123050612 tt a 2 + 25995895 a 3 - 133000 of - 2915745 a 2 /?)] . (59) 

We now suppose that the denominators in these three expressions for $n are differ- 
ent from zero. Later we consider the cases in which each of them is different from zero. We 
also suppose that spin coefficients a, (3, tt are different from zero. It was shown in CM that 
if one of these spin coefficients is equal to zero, the others must be zero too, and Huygens' 
principle is violated. 

The next step consists in proving that l|53|). and (|5*7|) imply that a, (3 and tt 

are proportional to each other. In order to get a system of with only two complex variables, 
instead of three, new variables are defined as follows: 

X\ := — , x 2 := — ■ (60) 

TT TT 

By subtracting from (|55[). taking the numerator and dividing by (8— X2—5x±) (5776^2 2 vr7f), 
we get 

N x := 178100 x x x 2 xf + 284960 x 2 xf + 208240 xf xl + 130150 xlxf x x 
+ 109825 xlxi x\ + 252850 x 2 xi x\ + 523744 x 2 xi + 3451480 x x xi 
+ 341900 xixi xl + 265888 xixl + 735800 x 2 x^xi + 2915264 x7 
+ 1018400 xix\ - 879008 xl + 18263488 - 408050 x\ xl + 4864450 x\ 
+ 35335248 x x + 22731900 x\ - 2101032 x x xl - 1622550 x\ xl = . (61) 

Subtracting (|55j) from (|57|). taking the numerator and dividing by (8 — x 2 — 5xi)(1527P7f), 
gives 

N 2 := -11651821200 x\ - 26531539120 x\ - 10132263424 - 26800626944a;i 
+ 242619584x2" + 593671488xl 2 + 2256829184x1- 35155250xJ 2 x 2 Xl 2 
+ 21550100 xf x 2 x x - 128589500x2 x 2 x\ + 11036720 xfxl 
+ 21755825 xf x\ xl + 130839250 x\ xi - 216631750 xix\ xl 

- 17700400 x^xi + 372958500 x\xl+ 220600700 xf x\ + 34480160 x 2 xT 3 

- 1915591500 x\ - 28320640 xf + 6897950 xf x x xl + 724005800 xf x x 

- 112640320xix 2 xi 2 - 90878112 x 2 xf + 59839936 xfxl 

+ 72209280 xTxT 2 x x - 1002142990 xlxixj - 585808600 x 2 xix\ 

- 448045312 x 2 xi + 632690016 x x xi- 1546809184 xix x xl - 796702240 xi xl (62) 

- 887976960 x 2 x2"xi + 509836460 xix\ + 2328634300 x\xl 

+ 5728848896 x x xl + 5470002280 x\ xl = . (63) 
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At this point we shall consider x±, X2, x±, ~x~2 as independent variables, and use the 
package grobner in Maple, for the polynomial system formed by the polynomials N±, N2. 
This package computes a collection of reduced (lexicographic) Grobner bases corresponding 
to a set of polynomials. The result is a list of reduced subsystems whose roots are those of 
the original system, but whose variables have been successively eliminated and separated 
as far as possible. In the present case we obtain four subsystems given by 

d := [42-3x1 + 65x2X2, 8 + 5xi ], 

G 2 := [205049562510x1x2 - 2072817918600 x 2 x^ + 529175067720 x 2 
+ 5001500073283 xf - 3239213905470x1+ 3029503111800 x^ 

- 26163100475032, 23707187714600 X2X2" 2 - 1207011 1345240 x 2 x2" 

- 15004500219849x1- 34648966659800 xf + 11975391986580x2" 

+ 41386627076564, -1113092 - 431311 xT+ 1909780x2"+ 954890x^x1, 
205xi + 368], 



G 3 := [2175607695654600868570x1x2 - 244429060944194171242925 x 2 X2"xi 

- 362016456337543432617920x2X2" - 25492004395136420363950 x 2 xi 

- 33777552239002428460240x2 - 352210319977170626297190 xixT 2 

- 527515033185400238012371 xf - 372609773697867989940085 xixT 

- 568758266358009596992694x1+ 357242473687668404124275 xix^ 
+ 529100974647178863056960 xi+ 421165196163010815629650 xi 

+ 613523694569903050334320, 74421671368200 x 2 xf 

- 372108356841000 x 2 x2"xi - 595373370945600 X2X2" 

+ 202992871981785xixI+ 309188233840256xT- 108770135076600 xf 
+ 408707348737250 x x xi + 731690446259960 x^ - 363026773505180 x x 

- 558465136160528, 

139740 X2" xT - 497365 xi xj - 799324 xT + 279480 xi - 1 187280 xi - 1879248, 
43975 x 2 + 137900xi + 107824] , 

G 4 := [15138500xTx 2 x2"- 200682625 x 2 x 2 xT- 589775940 x 2 xixT 

- 425769864 xlx 2 + 30277000 x 2 xi - 677199250 x 2 x\ - 2061417280 x 2 xi 

- 1553711328x2 - 34073270 x^xT 2 - 2202762525 x 2 xT 2 - 6761833290 xixT 2 

- 5193018500 xT 2 - 13933181776 - 6178746050 x\ - 18555219840 x x 

- 180544080x7- 17865753288x1- 158418580x7x1- 23445716600x1x1 

- 7696469075 x\ xT, 38482345375 x \ xi + 30893730250 x\ 

+ 11013812625 x\xf + 1003413125 x 2 x? xT + 3385996250 x 2 x\ 

+ 51211999525 x\ xf + 16053417900 x 2 x\ + 4651759450 x 2 xi 2 xl 

+ 177479507850 x\ xT + 140245780600 x\ + 212140488080 xi 

+ 7133355840 x 2 x x xT + 272934958520 x x xT + 25260582880 x 2 x 1 

+ 79415365840 xixT 2 + 106967929856 + 13183919424 x 2 + 41078949112 xf 
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+ 139995956352x1+ 3612843312x1x2, 427238747000 x 2 xf 

- 80593740760 x 2 - 961615825940 xf xl + 72041680708 x^xl 
+ 31587988349750 x\ + 9360849735875 x\ xl + 90434478667240 x Y 
+ 26404464403460xixT+ 549197297800xiX2"+ 64702426970096 
+ 18636073790528x1+ 1397883090496x2"- 2547657512880X2" 2 , 
31185310 x 2 X2" xi + 52924196 x 2 xi - 6814654 xi xT - 716849880 x\ 

- 440552505 x\ xl - 2156806048 x x - 1352366658 x 1 xj - 45578530 x 1 xi 

- 1620399064- 1038603700x1- 90980056x2", 505750 x^ + 149875 x\xl 
+ 1539520xi + 440460 xixT+ 116450 X2"xixT+ 232900xiX2"+ 1160352 
+ 317976x1+372640x2"+ 186320 xixl] . 

The only sets where the solutions x\ = const., X2 = const, are not obvious are G\ and 
G4. For G\, if we substitute x\ = —8/5 into iVi we get 195x2^2 + 702/5 = 0, which is 
incompatible with the first equation of this set. 

Let us consider now the fourth and fifth equations in set G4 given, respectively, by 

31185310 X2:x 2 xi + 52924196 X2"x 2 - 6814654 xj xl - 716849880 x 1 2 

- 440552505 x x 2 xl - 1352366658 xi xj - 2156806048 xi - 45578530 x x x$ 

- 1038603700x1 - 1620399064 - 90980056x1 = (64) 

and 

505750xi 2 + 149875xi 2 xT + 440460xixT+ 116450x2" xixl+ 1539520xi 
+ 232900xi X2"317976xT + 186320 x^ xl + 1160352 + 372640x2" = . (65) 

By applying grobner to (|64|) . (|65j) and their complex conjugates (in this case, two new 
equations), we obtain a system of polynomials for which all solutions have xi and X2 con- 
stant. 

Let us consider now the special cases in which the denominators in the previous 
expressions for 3>n are zero. The denominators of (|55|) . (|53|) . and (|57|) are given, respectively, 
by 

di := 10868X2" 2 + 772320 x^ + 1020276 x x xi + 4085 xf x x 

- 1158240 x\ - 2048352xi - 214700 x\ + 335985 xix\ - 1191680, (66) 

d 2 := 148960 + 42940 x\ + 10412 xf - 18564x2" - 13047 x x xi + 162944 x x , (67) 

d 3 := ( 148960 + 42940 x 2 + 10412 xf - 18564x2" - 13047 x x xi + 162944 x x ) 
(-24320 xf + 7684576x2" + 9460852 x x xi + 2915745 xi xf - 194465152 x x 

- 123050612 x\ - 102596352 - 25995895 xf + 133000 xf x x ) . (68) 

Applying 5 to (|6U)) we obtain 

63729588 xix\ - 670120 x 2 xa" 2 + 12371205 x\ xi - 617652 xf + 18435168 x\ 
+ 108263616xi xa"+ 10424160 X! 3 + 10725120 x x - 1048608 xf + 60746496x2" 
- 383325 xf xi - 1749232 xf x x + 1932300 x\ = . (69) 
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(70) 
(71) 



Applying 5 again on (|69|). gives 

- 71849032 xf x x - 2917197592 xf xf + 2126429184 xi x\ 
+ 3683108352 xi xi + 357185400 x\ xi- 128701440 x\ - 125089920x1' 

- 67408896 xi 3 - 18590290 x\ xf + 6696360 xi 4 + 4179810 x x xf 

- 221222016 x\ - 2418547200 xf + 2059223040x1 - 22411650 xf xi 

- 614629870 x\ xf - 23187600 x\ - 4605932928 xf x x = . 

Using grobner package on (|66jl. (Jfi9"|) and (|7T7|) gives the empty set solution. Applying 5 to 
(joT|) . gives 

- 128832 xf - 101344 xf x x - 2591852 xixf - 8248048 x x xi - 257640 xf 

- 6550592 xi - 977664 xf - 20824 xf - 893760 x a = . (72) 

Applying 5 again on (J72J), gives 

62472 xf - 60094064 xi x\ + 35714228 xf xf + 2832764 xf xi 

- 2838720 xf + 8043840 xf - 210698752 x x xi + 8798976 xf + 86775744 xf 

- 171601920 xi - 1690904 xf x x + 111204048 xf x x + 2318760 x\ = . (73) 

Using grobner package on (|67|). (f?2]) and (f?3|) we obtain the empty set solution. We 
observe now that one of the factors in dj, is di. Thus, if ds = 0, we need to consider just 
the expression 

133000 xf xx - 24320 xf + 7684576 xi + 9460852 x x xi + 2915745 xi xf 

- 194465152 x x - 123050612 xf - 102596352 - 25995895 x\ = . (74) 

Applying <5 twice on this equation, results in 

1923742456 xi xf - 54838615 xf xf + 387128860 xf xi + 28673280 xi 3 
+ 1750186368 xf + 3181762656 xi xi + 1107455508 xf + 923367168 x x 

- 143083296 x^ 2 + 1751900928 xi + 17772600 xi 3 xi - 175990444 xi 2 x x 

+ 233963055 xf = (75) 

and 

3071183072 xf x x - 65381238064 xf xi 2 + 363466389312 xi xf 
+ 414477092352 x x xi + 141775882688 xf xi - 11080406016 xf 

- 13289466096 xf + 2839158528 xi 3 + 811332320 xf xf - 315187200 x 2 4 

- 196695600xi xf - 21002236416 xf - 58327724544 xi 2 + 177286496256 xi 
+ 20770389380 xf xi - 13320040240 xf xi 2 - 2807556660 xf 

- 106950649152 xi 2 x x = . (76) 
Using grobner package to (f73|) and (|76*|) we obtain again the empty solution set. 
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We need now to study the case in which x\ and X2 are constants. From Sx% = 
d(a/w) = and 8x2 = 5(f3/w) = we get, respectively, 

31xi + 10 xi 2 + 2 = (77) 

and 

63x 2 + 40x 2 xT + xT = 0. (78) 

The above equations have two solutions, given by x\ = —3/2, x 2 = 1/2 and x\ = x 2 = —8/5. 
The first one satisfies 5xi + X2 — 8 = 0, which will be considered next. The second case is 
impossible, since these values don't make N\ equal to zero. 
Let us now consider the case: 



- 5a - 8vr = , (79) 
or, using variables xi and x 2 , and dividing by ir, 

8 - x 2 - 5xi = . (80) 
Applying 5 to this equation, using (|3*2*)l . (|STj) and 1)52(1 . we get 

34x 2 xl - X2 + 56x 2 + 15xT 2 + 24x1 = . (81) 
The only solution for both equations is given by 

xi = -3/2, x 2 = l/2. (82) 

Since the numerator on the right side of (|48|) must be zero, we obtain, solving for <J>n, 

86 

$11 = "141^- (83) 
Applying 5 on this equation, we obtain 

Sir = -7r?f . (84) 

Equations (JSHJ) and (JHU give 

5tt = tt 2 . (85) 

Computation of the commutator [5, 5]ir, using (|84*|) and (|85j) gives 7r = 0. 

Finally, we consider the case $n = 0. Here both denominators of (|55|) . and of (|57(1 
must be equal to zero. Applying Buchberger-Grobner algorithm to them and their complex 
conjugates and using the variables defined in (|60|) . we can verify that the only possible 
solutions again fall in cases we have studied before, i.e., either xi or X2 are zero, or both xi 
and X2 are constants. 
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5 Conclusions 



Theorem 1 was proved for the case of Maxwell equations, i.e., there are no Petrov type III 
space-times on which Maxwell's equations satisfy Huy 'gens' principle. For the neutrino case, 
k\ = 8, £?2 = 13, the proof is similar |18| . The use of the package NPpsinor in Maple was 
essential for the conversion of the Alvarez- Wiinsch five-index necessary condition from the 
tensorial to dyad form. 

The polynomial system obtained from the integrability conditions was simplified 
using Maple's package grobner. Since a direct application of the algorithm is apparently 
impossible due to the large size of the polynomial system, a "divide and conquer" approach 
was applied with success, showing that the necessary conditions III and V for the validity 
of Huygens' principle cannot be simultaneously satisfied for Maxwell's equations in Petrov 
type III space-times. 
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